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Abstract 

We study the ellipticity and the "Nekhoroshev stability" (stability properties for finite, but very 
long, time scales) of the Riemann ellipsoids. We provide numerical evidence that the regions of 
ellipticity of the ellipsoids of types II and III are larger than those found by Chandrasekhar in 
the 60's and that all Riemann ellipsoids, except a finite number of codimension one subfamilies, 
are Nekhoroshev-stable. We base our analysis on a Hamiltonian formulation of the problem on 
a covering space, using recent results from Hamiltonian perturbation theory. 



1 Introduction 

A. The Dirichlet problem consists of the "linear" motions of an ideal, incompressible, homogeneous, 
self-gravitating fluid mass having an ellipsoidal shape at any instant. This problem originated in 
attempts to determine the shapes of the planets, and has a long and distinguished history (see 
e.g. chapter 1 of || and references therein). Let us only remark here that Newton and Maclaurin 
established the existence of rigidly rotating oblate spheroids and that Jacobi described rigidly 
rotating asymmetric ellipsoids. Dirichlet proved the existence of "linear" fluid motions, such 
that the position x(t, y) at time t of any material particle y is given by 

x(t,y) = F(t)y, F(t) €SL(3). (1.1) 

If, as we assume, the reference configuration of the fluid mass is a ball of radius p, then the free 
surface of the fluid determined by the motion (|1.1| ) is at every instant an ellipsoid with semiaxes 
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pai, pa>2, pa^, where a\, 0,2 and 03 are the singular values of the matrix F (namely, the eigenvalues 
of VfW). 

The Maclaurin spheroids and Jacobi ellipsoids are steady rigid motions of the Dirichlet problem. 
Dedekind identified a new class of steady, but non-rigid, motions of the problem; they were con- 
structed by transposition of the Jacobi ellipsoids, thus bringing the particle relabelling symmetry 
of the system into play. Riemann reformulated the equations of motion of the Dirichlet problem 
in a particularly convenient form for studying asymmetric configurations (as we shall show, this 
reformulation corresponds to the passage to a four-to-one covering space). Riemann also classified 
all steady motions of the system, showing that they are generated by combinations of spatial and 
body symmetries, and studied their Lyapunov (or nonlinear) orbital stability. 

The steady motions found by Riemann, which are called Riemann ellipsoids, are motions ( |1.1| ) 
such that 

F(t) = exp(tf^) A exp(-tft r ) 

for some constant diagonal matrix A = diag (ai, 02, 03) and constant antisymmetric matrices 0/ 
and £l r . Except for the special case of the Jacobi ellipsoids, for which D, r = 0, Riemann ellipsoids do 
not perform rigid motions: the motion is the composition of an internal rotation, a stretch along the 
principal axes, and a spatial rotation; the free surface retains a rotating ellipsoidal shape, but the 
fluid particles describe rosette-shaped motions which are either closed (periodic) or open (quasi- 
periodic) depending on the two angular frequencies oj\ and uj t (that is, the vectors corresponding 
to the antisymmetric matrices Qi and Q r ). Riemann proved that there are five types of Riemann 
ellipsoids with distinct semiaxes, characterized by different relations between the angular velocities 
and the semiaxes. For two types (hereafter referred to as S2 and S3), both uji and u> r are parallel to 
the same principal axis of the ellipsoid; this axis is either the shortest (S3) or the middle (S2) one. 
For the other three types, called I, II, and III, both ui and iv r belong to one of the two principal 
planes containing the longest axis of the ellipsoid]]] 

Since the Riemann ellipsoids are reduced equilibria, namely equilibria modulo symmetries, of 
the Dirichlet problem, it is meaningful to study their stability modulo symmetries, which we will 
frequently refer to simply as stability. Riemann found that all ellipsoids of type S3 and a subfamily 
of the S2 ellipsoids are stable. Riemann's analysis does not lead to any definite conclusion about the 
(in)stability of the remaining ellipsoids, because they are saddle points of his Lyapunov function. 

More recently, Chandrasekhar studied the 'ellipticity', or 'spectral stability', of the Riemann 
ellipsoids. We say here that a (relative) equilibrium is elliptic if the eigenvalues of the linearization 
at the equilibrium are purely imaginary; note that, accordingly, we treat the case of zero eigenvalues 
as elliptic, even though this is not customary. Chandrasekhar [UlUH numerically found a region 
of ellipticity of the ellipsoids of type S2 larger than the region of Lyapunov stability found by 
Riemann. He also found nonempty regions of elliptic ellipsoids of type I and III, but none of type 
II. As we shall see, some of these conclusions are not entirely correct, and in fact require substantial 
modification. 

B. For a Hamiltonian system such as the Dirichlet problem, ellipticity is a necessary, but not 
sufficient, condition for Lyapunov stability. Therefore, the ellipticity analysis is of twofold interest: 
on one hand, the absence of ellipticity implies Lyapunov instability; on the other hand, it opens 
up the problem of the stability properties of those Riemann ellipsoids that are elliptic, but of 

1 The distinction between S2 and S3 ellipsoids that we make here is not standard; Chandrasekhar calls them 
collectively ellipsoids of type S. 
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undetermined Lyapunov stability. The present article is mainly devoted to a study of the latter 
question. 

First of all, as a prerequisite for this analysis, we repeat the Lyapunov stability and the ellipticity 
analysis of the Riemann ellipsoids. While our analysis confirms Riemann's Lyapunov stability 
results, we find significant discrepancies with Chandrasekhar's results: there is a nonempty region 
of elliptic ellipsoids of type II and a region of elliptic ellipsoids of type III that is substantially larger 
than the one reported by Chandrasekhar. Therefore, the regions of known (Lyapunov) instability 
of the ellipsoids of types II and III are substantially smaller than those found by Chandrasekhar. 
Moreover, the ellipticity regions of the ellipsoids of type I, II and III have a finer and more complex 
structure than was previously realized. 

The main goal of the article is to use the methods of Hamiltonian perturbation theory to 
investigate the long-time stability properties of the elliptic Riemann ellipsoids of undetermined 
Lyapunov stability. In fact, the Dirichlet problem is a Hamiltonian system with symmetry and 
the Riemann ellipsoids are the relative equilibria of this system. They correspond to the equilibria 
of the reduced system, which has four degrees of freedom for generic values of the momenta and 
three degrees of freedom in special cases. (The reduced phase space is in fact an orbifold, with a 
singular set containing the Riemann ellipsoids of type S, but this difficulty is overcome by passing 
to Riemann's covering space.) Since the reduced system has more than two degrees of freedom, 
KAM theory cannot be used to show the stability of its equilibria. Therefore, we shall base 
our investigation on Nekhoroshev theory, which allows one to obtain "practical" stability results, 
namely, to prove that motions starting sufficiently near the equilibrium remain near it for times 
growing exponentially with (some power) of the initial distance; specifically, we shall use some 
recent results obtained in (see also 17], |ll|). The results of our investigation indicate 

that all elliptic Riemann ellipsoids are Nekhoroshev-stable, with the possible exception of a finite 
number of codimension one families of ellipsoids satisfying certain low-order resonance conditions. 

In view of the complexity of the system, we resort to numerical methods at certain steps in 
the ellipticity analysis and in the construction of the normal forms for the reduced Hamiltonian. 
We numerically compute the eigenvalues of the linearization and numerically construct the normal 
forms for a finite (though rather large and, we think, significant) number of values of the momentum. 
Therefore, the forementioned results are not rigorous, but we think that our analysis provides strong 
evidence for our conclusions. 



C. Our analysis is based on a Hamiltonian formulation (on a covering manifold) of the Dirichlet 
problem, which is briefly developed in Section 2. For the sake of brevity, we skip most of the 
computational details and focus on the conceptual features. The Riemann ellipsoids are described 
in Section 3. Overall, Sections 2 and 3 contain a treatment of the Dirichlet problem which, although 
concise, is complete and detailed; this seems necessary, because the existing literature does not seem 
to include a comparably comprehensive treatment. 

Sections 4-6 contain the ellipticity and the Nekhoroshev-stability analyses. In Section 4 we 
review the results on Lyapunov stability and instability and describe the results of our ellipticity 
analysis; a detailed comparison with Chandrasekhar's results is included in Appendix B. We briefly 
review the necessary notions for the Nekhoroshev-stability analysis in Section 5 and describe the 
results for the Riemann ellipsoids in Section 6. A Conclusion follows. 

All of the numerical computations have been performed using the software package Mathemat- 

ica. 
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2 The Dirichlet Problem 

A. Dirichlet's equations. Dirichlet § showed that the linear motions fll.l| ) form an invariant 
submanifold of the phase space of an ideal, incompressible, homogeneous, self-gravitating fluid, 
with the boundary condition that the pressure is constant at the free surface. He expressed the 
restriction of the hydrodynamic equations of motion to this submanifold as equations for the curve 
F(t) E SL(3), obtaining a second order system on SL(3). 

To describe these results, we introduce the following notation, which will be used throughout 
the article. We denote by • and || || the Euclidean scalar product and norm on R 3 and by D p the 
ball {y E M 3 : \\y\\ < p}. We use the standard inner product and norm on L(3) given by 

(A, B) = tr(AB T ) and \A\ 2 = (A, A). (2.1) 

(Note that the subspaces of symmetric and antisymmetric matrices are ( , ) -orthogonal comple- 
ments of one another in L(3).) We shall frequently identify M 3 and so(3) by means of the isomor- 
phism w ^ (i, where u>£ = to x £ for all £ E M 3 . Note that (u>,£) = 2uj ■ £ for all u, £ E M 3 . 
Following Riemann, we utilize the "singular value decomposition" of matrices, nowadays a standard 
tool in numerical linear algebra: If F E SL(3) and oi, 02, a^ are its singular values, then there exist 
matrices Ui,U r E SO (3) such that 

F = UiAUj , A = diag (ai, o 2 , a 3 ) . 

We shall refer to the triplet (A, Ui, U r ) as a "SVD" of the matrix F. The matrices A, XJ\ and U r are 
not unique, but one can clearly fix A by ordering its entries, e.g. a\ > 02 > 03 (see also Proposition 
H below). Finally, we set 



r°° r 

V(a) = -2ng / {s + a?)(s + a 2 2 )(s + a\) 
Jo L 



ds , 



where g denotes the gravitational constant; by a slight abuse of notation, we shall regard V and its 
derivatives as functions either of a = (01, 02, 03) or of the diagonal matrix A = diag (ai, 02, 03). 

Dirichlet's result can now be expressed as follows: Consider a curve t F(t) C SL(3) defined 
on some open interval 3 C M and let (A(t),U[(t),U r (t)) be any singular value decomposition of 
F(t). Then the map 

x : 3 x D p — ► M 3 , (t,y) » F(t)y 

is a solution of the hydrodynamical equation for an ideal, incompressible, homogeneous, self- 
gravitating fluid with constant pressure at the boundary iff 

F + UiV'(A)Ur~^ = , (2.2) 

where V = diag |J, and ¥ F G = G-\{G, F)F~ T for any G E L(3). Note that the 

matrix U{V \A)Uj is independent of the particular SVD of F used; this is proven by observing that, 
ifU,Ve SO(3) and A = UA'V with A and A' diagonal, then V{A) = V'QJA'V) = UV'(A')V. 
In fact, Dirichlet used the polar decomposition, not the SVD, when writing the force term in the 
equations of motion. 

Remark: In his work on Dirichlet problem, Chandrasekhar required that the relative pressure 
p — p , where p Q is the constant pressure on the surface of the fluid, be nonnegative at any point 
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of the fluid. Following Riemann, we do not make this assumption. (There are many fluids that 
support even very high tensile states.) 

Since the gravitational forces are conservative and invariant under spatial rotations and particle 
relabelling, equations ( |2.2| ) inherit the conservation of energy, angular momentum and circulation 

2 

from the hydrodynamic equations. For future reference, we observe that if we set c = where 
m denotes the mass of the fluid, then the energy of any motion (|1.1|) is c[^|F| 2 + V(a)] and the 
angular momentum is the vector k satisfying k = 2c skew(FF T ); here skew(G) = ^(G — G T ). 
The conservation of circulation is equivalent to the constancy of skew(F T F). 



B. Riemann's equation. Shortly after Dirichlet's work appeared, Riemann [20] rewrote equation 
(|2.2| ) using what is nowadays called the singular value decomposition of matrices and left-trivializing 
the factors TSO(3), namely expressing the tangent vectors Ui and U r in terms of the antisymmetric 
matrices = Uf(j l and £l r = U^U r . Once these substitutions have been made, equation 
takes the form 



A + 2(QiA - An r ) + (n t A - An r ) + (nfA-2n l An r + Anf) + v'(A) 







(2.3) 



(which should be compared with equations (a) in pOfl ), This equation determines a second order 
system on A x SO(3) x SO(3), where 



A 



diag (oi, a,2, (aia 2 ) 



: a\ > a,2 > a 



(if A has repeated eigenvalues, 0/ and Q r are not uniquely determined by (|2.3j )). This system is 
equivalent to the restriction of Dirichlet's equation (p.2[) to the submanifold 



Q = {F G SL(3) : F has distinct singular values} . 

More precisely, because of the nonuniqueness of the SVD, the two systems are related not by a 
diffeomorphism, but by a four-to-one covering map: 



Proposition 1 The map 

e:Ax SO(3) x SO(3) -»■ Q , {A,U h U r ) ^UiA 
is a ^-fold covering map which maps solutions of equation into solutions of equation ([2.fy . 

Proof. If (A,Ui,U r ) and (A,Vi,V r ) are any two SVDs of F G Q, with A = diag (cti, «2, 0,3), 
a\ > a.2 > 03, then V\ = B4U1 and V r = RiU T for some i = 0, . . . , 3, where 

Ro = l, R x = diag (1,-1,-1), R 2 = diag (-1,1,-1), and R 3 = diag (-1, -1, 1). (2.4) 

Thus, the map C is four-to-one. That it is a submersion is proven, for instance, by showing that 
the determinant of its Jacobian is nonzero. (Using 01,02 as coordinates on A and left-trivializing 
the factors TSO(3), one obtains the determinant {a\ — a^X ? — °i)( a 2 — a i)-) B 

In Riemann's formulation, the use of the SVD has the advantage that the absence of the matrices 
U[ and U r from the equation ( |2.3j ) makes manifest the SO(3) x SO(3) invariance of the Dirichlet 
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problem. Prom our perspective, the formulation on the covering also has the advantage that, once 
we pass to the Hamiltonian formulation, the symplectically reduced phase space associated to the 
cotangent bundle of the covering is a manifold, while reduction of the system on T*Q produces 
singularities. 

Remarks: (i) The submanifold Q is not invariant under the dynamics of the Dirichlet problem 
(axisymmetric configurations can evolve into asymmetric ones, and vice versa) but this is not a 
problem for our analysis: the Riemann ellipsoids with asymmetric configurations clearly belong to 
Q for all time and, since Q is open, it will be a consequence of our stability estimates that nearby 
motions also remain in Q for the time under consideration. 

(ii) The restriction to the case of distinct singular values excludes all axisymmetric Riemann 
ellipsoids. The stability properties of the Maclaurin spheroids, which rotate rigidly about the axis 



of symmetry, are completely known (see [20|||[14]); the remaining axisymmetric ellipsoids merit 
further investigation. 



C. Hamiltonian formulation. We now describe the Hamiltonian formulation of Riemann's 
equation ( |2.3D on the cotangent bundle of the covering manifold A x SO(3) x SO(3). We shall use 
the first two diagonal entries a\ and a2 as coordinates on the set A, by means of the diffeomorphism 
a : 23 — ► A, where 

£ = {(hM) G ■ h > b 2 > &7 1/2 } and a(h, b 2 ) := diag (h, h, ^ 

(We denote by b = (61,62) the first two singular values, so as to distinguish b G K 2 from a£l 3 ; 
wherever we write 63, it stands for (fei^) 1 , riot for an independent variable.) Furthermore we 
identify the cotangent bundle of each factor SO(3) with SO(3) x R 3 as follows: we first identify the 
cotangent space Tj}SO(3) with T^SO(3) by means of the inner product |( , ) (see (2.1)) and then 



identify Tj/SO(3) with ~R S by means of the " isomorphism, thus associating to each jj, G T^SO(3) 
the unique m G M 3 satisfying 

1 3 
fi-Uu = —(m,Lu) = m-oj V uj G R . 

After these identifications, we work on the sixteen-dimensional manifold 

M = Bxl 2 x (SO(3)) 2 x (R 3 ) 2 , 

which is diffeomorphic to T*(A x SO(3) x SO(3)). We shall denote by (b, c, U,m) the elements of 
M, with b G S, c e M 2 , U = {U h U r ) G (SO(3) 2 ), and m = {m u m r ) G (M 3 ) 2 and use the induced 
scalar and vector products on (M 3 ) 2 given by 

(u h U r ) ■ (&,f r ) = Urtl + Wr-Zr, {^h ^r) X tr) = fal X Zh u r X tr) ■ 

When writing tangent vectors to M we shall again left-trivialize and identify the tangent spaces 
to the two factors SO(3) with IR 3 ; thus, (v,w,ui,n) denotes the vector (v,w,(UiUi,U r u r ),n) G 
T(b,c,u,m)^" Finally, in order to avoid the introduction of a new symbol, we write V(b) to denote 
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Proposition 2 Riemann's equation (2.S) on A x S0(3) x S0(3) is equivalent to the Hamiltonian 
system on the manifold M defined by the Hamiltonian 



H(b,c,U,m) = 4c-K(6)c + \ m ■ J(b)m + V(6) , 



where 



K(6) 
J(6) 



1 



6? (&I + « 



/J! (6) J 2 (&)\ 
U(&) Ji(6)J ^ 



-6 3 
&| (6? + 61) 



J! (6) - diag (^|j2 , (^2)2 

TfM - Hi.no- ( gMg 2fe i fe 3 2bi6 ? 



(2.5) 



and i/ie "left-trivialized" (in the sense just specified) symplectic structure a is defined by 
a(b, c, U, m)((v,w,uj,n), (v ,w',cj',n )) = v-w 1 — v ■ w + u ■ n — n-u + m-(uxw). (2-6) 



This proposition can be proven by a direct computation showing that Hamilton's equations 
for the Hamiltonian system defined by ( |2.5| ) and (|2.6| ) are equivalent to Riemann's equations. 
Alternatively, one can observe that Riemann's equations are equivalent to the Lagrangian system 
defined on TA x TS0(3) x TS0(3) by the Lagrangian 

H ((A, A), Ui&i, U r oj r ) = ±|i + CjiA - A£o r \ 2 - V(A) (2.7) 

and construct the Hamiltonian system using the Legendre transform. For the sake of brevity we 
skip all details, remarking only that the relation between 'angular momenta' and 'angular velocities' 
is uj = J{b)m. 

Remark: The Lagrangian nature of Dirichlet's problem, with Lagrangian equal to the difference 
of the kinetic and the potential energies, is implicit in Riemann's work and has been repeatedly 
used since then. One can verify by a direct computation that the Euler-Lagrange equations for the 
Lagrangian ( |2.7| ) are equivalent to Riemann's equations (2.3). 



The Hamiltonian system (M, H, a) of Proposition ^ is invariant under the symplectic action $ 
of S0(3) x S0(3) on M given by 

®(R u R r ){b,c,{U u U r ),m) = (b,c,(RiUi,RrU r ),m). 

Since the group does not act on the factor 23 x M. 2 , while it acts on each factor S0(3) x M 3 as the 
left-trivialized cotangent lift of the left action of S0(3) on itself, all results from the elementary 
case of the left action of S0(3) on T*S0(3) apply (see, e.g., or |?J). Thus, one concludes that 
the action $ has the momentum map 

J : M -> (M 3 ) 2 , (a, c, U, m) (^mj , U r m r ) . 

For any r\ G (M 3 ) 2 , one may identify J -1 (77) and 23 x K 2 x (SO(3)) 2 , in which case the immersion 
i Tj : J _1 (?7) —>■ M is given by i v (a,c,U) = (a,c,U,(jffr]i,Ujr] r )). The two components Uimi and 
U r m r of the momentum map coincide, up to constant rescalings, with the total angular momentum 
and the circulation; see Section 2. A. 
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The next Proposition, which describes the reduced system (P„, a^, H^), also follows immediately 
from standard results for the reduction of an SO(3) invariant Hamiltonian system on T*SO(3). We 
denote by the sphere of radius p. As usual, G v denotes the isotropy subgroup of r\ with respect 
to the coadjoint action; in the case at hand, G v = {rotations about r]i} x {rotations about r/ r }. 



Proposition 3 (i) For any r) = (r)i,7]r) £ (M 3 ) 2 with iji / / rj r , the reduced phase space 
J~ 1 (r])/G ri can be identified with 

P v = Sxl 2 x(4nx4 r|| ). 

After this identification, the canonical projection tt ti : J _1 (r/) — ► P v is n v (b, c, U) = (6, c, (Ufrji, Ujrj r )) 
The symplectic structure on defined by the equality z*cr = ir^a^ is 

CTn(b, c, m)((v, w, m x u), (v , to, m x &)) = v ■ w — v ' • w — m • to x to' . 

The reduced Hamiltonian H v : P v — > M, defined by H v o 7r^ = H o i v> is 

Hr,(b,c,m) = ic-K(6)c + \m ■ J(6)m + V(b) . (2.8) 

(ii) If rj = (rji,0) or r] = (0,r]i) with rji ^ 0, then the reduced phase space can be identified with 
Bxl 2 x S| .I and one has ir, n (b,c,Ui) = (b,c,Ufrn) and 

a v (b,c,mi)((v,w,mi x lui), (v' ,w' ,mi x u![)) = v ■ w' — v' • w — mi • oj\ X u[ 
H v (b,c,mi) = ic-K(b)c + \mi-J x (b)mi + V(6) . 

The reduced system generically has four degrees of freedom, but only three in the special case 
(ii), which arises in the case of the so-called "irrotational" Riemann ellipsoids. The case rji = rj r = 
yields a reduced system with two degrees of freedom, the only equilibrium of which is the stationary 
sphere. 

In view of the application of Nekhoroshev theory, we note that the reduced Hamiltonian H v is 
an analytic function. This is seen by observing that the self-gravitational potential can be written 
in the form 

where F(ip \ k) is the incomplete elliptic integral of the first kind, which is analytic for ip G (0, ^) 
and k G (0, 1) (see @, pg. 4-5, 299). 

D. Equivalence and symmetry of reduced systems. As we now discuss, (i) certain classes 
of reduced systems are equivalent and (ii) the reduced systems have certain discrete symmetries. 
For the sake of conciseness, we consider here only the generic case, in which both r]i and rj r are 
nonzero. (The case in which one of the momenta vanishes is recovered with obvious modifications.) 
There are two distinct reasons for equivalence of different reduced systems: 
El. The reduced system (P^, H v , a v ) depends on -q = (r]i,rj r ) only through the norms of r\i and rj r . 
This equivalence reflects the invariance of the Dirichlet problem under spatial rotation. 
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E2. Dedekind proved that if F(t) is a solution of the Dirichlet problem, then F(t) T is also a solu- 
tion [||. (The proof follows immediately from the observations that ¥ f tG t = (¥pG) T for any 
matrix G and that \J\AUj is a singular value decomposition of F iff U r AUf is a singular value 
decomposition of F T .) Chandrasekhar calls 'adjoint' the two motions determined by F(t) and 
F{t) T . This invariance of the Dirichlet problem has a counterpart in the Hamiltonian formula- 
tion on the covering space M, which is invariant under the symplectomorphism exchanging XJ\ 
with U r and mi with m r . This map takes solutions to solutions, but exchanges the level sets 
of J. Correspondingly, the reduced systems for the values (rji, rj r ) and (r/ r , rji) of the momentum 
map are conjugated to one another by the diffeomorphism (b, c, (mi,m r )) \— > (b, c, (m r ,mi)). 
Furthermore, each reduced system is invariant under two discrete actions: 

i. The Z2~action (b,c,m) t— > (6, c, — m), which reflects the invariance of the original system with 
respect to the combination of spatial reflection and time reversal. 

ii. The Z4-action (b,c,m) i— ► (b,c, Rim), i = 0, ...,3, where Ro,...,R^ are given by and 
Rj(mi,m r ) = (Rim i, Rim r ). The points (b,c,Rim) are obtained from one another by a re- 
flection about one and the same coordinate axis on the two spheres; see figure [|. Hence the 
Z4~orbit of the point (b, c, m) consists of two distinct points if m is parallel to a basis vector 
(ej, ej), j = 1, 2, 3, and of four distinct points otherwise. 

In the stability analysis, we need not distinguish equivalent reduced systems, nor distinct points 
within the Z4 x Z2 orbits, in each reduced system. A Z4 x Z2~orbit consists of eight, four, or two 
points, depending on the number of nonzero components of mi and m r . We will say that all points 
in such an orbit are Z4 x Z2-equivalent and regard two points in the reduced space as adjoint if 
any two points in their Z4 x Z2~orbits are adjoint. 

By a Riemann ellipsoid we will mean an equivalence class, under spatial rotations (El) and 
transposition (E2), of motions of the Dirichlet problem corresponding to an equilibrium of the 
reduced system. Thus each Riemann ellipsoid can be identified with the Z4 x Z2~orbit of an 
equilibrium of some reduced system .0 

Remark: The Z4-symmetry is not present in Dirichlet's system Q2.2|) , but is clearly inherited 
from the passage to the fourfold covering manifold M of T*Q. (The existence of this symmetry of 
Riemann's equation is mentioned by Chandrasekhar (j|], pg. 72), but its origin is not identified 
there.) Since all points in a Z4-orbit correspond to the same state of the fluid, the passage to the 
covering introduces some redundancy, but allows us to avoid singular reduction. In fact, Dirichlet's 
equations ( |2.2[ ) determine a SO(3) x SO(3)-invariant Hamiltonian system on T*Q; however, the 
action is not free and the reduced phase space is the orbifold Sxl 2 x (S 2 x 5 2 )/Z 4 , rather than 
a smooth manifold. The singular set consists of all states in which mi and m r are both parallel to 
the same axis of the ellipsoid, and thus contains the Riemann ellipsoids of type S. 

3 Riemann Ellipsoids 

We now study the equilibria of the reduced systems of Proposition ||, and thus the Riemann 
ellipsoids. The results of this analysis are due to Riemann. As we mentioned in the Introduction, 
there are five classes of Riemann ellipsoids and the momenta mi and m r of the ellipsoids of each type 
are completely determined by their semiaxes. Correspondingly, in the set "B there are "existence 
regions" of the Riemann ellipsoids of each type. In order to specify these regions, and the momenta 

2 Note that, as was previously stated, we exclude from our consideration all axisymmetric ellipsoids. 
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Figure 1: Projection of a Z4-orbit on one of the spheres. 



of the corresponding Riemann ellipsoids, we introduce some notation. First we set 

-3/2 



C n (x,y,z) = 2ixg 
The functions 

G b ± (x,y,z) = 



[s + x z )(s + y z )(s + z z ) 



s n ds. 



n = 0,1,2. 



xz 



[[xy 2 z ± (x 2 y 2 - x 2 z 2 + y 2 z 2 )]C 1 (x,y, z) + (xz ±y 2 )C 2 (x,y, z) 



determine the momenta of the ellipsoids of type S, while the functions 

G(x,y,z) = x 2 (y 2 -z 2 )C x {x,y,z) + (y 2 - 4z 2 )[z 2 d(x, y, z) + C 2 (x, y, z)\ 



D(x, y, z) = x 2 (y 2 - z z ) + z z (Az 



G^(x,y,z) = {yTz) 



(y ± 2z) 



x 2 — z 2 G(x, y, z) 
x 2 — y 2 D(x, y, z) 



determine the momenta of the ellipsoids of type I, II and III; we refer to the latter three classes of 
ellipsoids collectively as type R ellipsoids] 3 ] Using these functions we define the five sets 



(3.1) 



and five pairs of real- valued maps one for each type a of ellipsoids; these maps and their 
domains are described in Table 1, where {e±, e 2 , e 3 } denotes the standard basis vector of M 3 and 

N±(x,y,z) = ■, 



Ss 2 


= {be-B 


G s (b l ,b 2 ,b 3 )>0} 






= {6e£ 


G s + (h,b 3 ,b 2 )>0} 




Si 


= {be £ 


bi < 2b 2 - b 3 } 




Bn 


= 


h > 2b 2 + b 3 , D(h,b 3 ,b 2 ) 


< 0} 


■Bin 


= {be £ 


h > b 2 + 2b 3 , G(h,b 2 ,b 3 ) 


>0} 



GUx,y,z) ± JGL(x,y,z 



t = R,S. 



3 In Riemann's article, the study of the existence of these ellipsoids precedes that of the ellipsoids of type S. 
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Type (a) 


Domain 




s 2 


3s 2 


A r ±(^i,&2,fc3)e 2 


S3 




W|(&i,&3,62)e 3 


I 


3i 


^(6i,63,62)ei + JV£(63,6i,62)e 3 


II 


3 n 


A^(&i,6 3 ,6 2 )ei+iV«(&3,6i,6 2 )e3 


III 


3m 


N£(h,b2, h)ei + N«{b 2 , h,h)e 2 



Table 1: The momenta of the Riemann Ellipsoids 



Proposition 4 (i) A point (b,c,m) £ P v is an equilibrium of the reduced system (P„, H n , a^), 
rjl 7^ 7^ r/ r , iff it is Z4 x ^-equivalent to either the point (6, 0, (/i+(&), or f/ie poini 
(6, 0, (fj,~(b), /i+(6))) /or some 6 G 23 a and some a = S 2 , S3, 1, II, III. 

(ii) If one component of r\ equals zero, then the equilibria of the six-dimensional reduced system 
are the points Z4 x ^-equivalent to (6, 0, fj,^(b)). This situation is met only at the equilibria of type 
S 2 satisfying G+(6i, 6 2 , 63) = G_(&x, 6 2 , 63) and those of type I satisfying b\b\ + b\b\ + b\b\ — 3b\ = 0; 
in &o£/i cases, n^i ) 7^ = Ma (^)- 

The proof of this proposition can be reconstructed from Riemann's and Chandrasekhar's work. 
A sketch of the argument is given in Appendix A. 



The existence regions (|3.1|) are drawn in figures || using the coordinates (x,y) = (o 2 /6i, 63/61), 
which have the advantage that the unbounded region "B is mapped onto the bounded triangle 
{(x, y) : < y < x < 1}. The borders of the various regions are determined by the definitions ( 0.1]) . 
Note that, since the existence regions of ellipsoids of different types have nonempty intersections, 
there exist Riemann ellipsoids of up to four different types having equal semiaxes. 

The Riemann ellipsoids belonging to a six-dimensional reduced space have either zero angular 
momentum or zero circulation and are called "irrotational" by Chandrasekhar. The existence of 
irrotational ellipsoids of type S 2 was known to Riemann, while the existence of those of type I seems 
not to have been noticed before. Note that the curve of irrotational ellipsoids of type S 2 divides 
the region 2>s 2 into two parts, in one of which the two momenta are coparallel, while in the other 
they are counterparallel.0 

Remark: Since Chandrasekhar requires that the relative pressure be positive inside the el- 
lipsoids, his region of existence of the ellipsoids of type II is smaller than ours, lying below the 
"zero-pressure curve" D (b\, 63, 62) Co + §b\C\ + 3C 2 = shown in figure ||.c (see ||, chapter 7, 
formula (195)). 



4 Stability and Ellipticity of the Riemann Ellipsoids 

A. Results. We now review the known results about the Lyapunov stability of the Riemann 
ellipsoids, regarded here as equilibria in the reduced phase space, and we describe the ellipticity 
analysis. Prom now on, we restrict ourselves to the generic case in which the reduced system has 
four degrees of freedom; the irrotational Riemann ellipsoids will be studied elsewhere. 



4 We say that two vectors u and v are coparallel (counterparallel) if u ■ v = ||w|| \\v\\ (= — ||m|| 
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t>3 :bi b3 :bi 




(c) (d) 

Figure 2: The existence regions of the five types of Riemann ellipsoids. The dashed curves in regions !Bs 2 
and ®i identify the irrotational ellipsoids and their adjoints (for which v a = 0); the former separates the 
coparallel (lower) and the counterparallel (upper) S2-ellipsoids. The dashed curve in region H>u identifies 
the zero-pressure ellipsoids. 




Figure 3: Regions of ellipticity of the Riemann ellipsoids of types I, II and III (light gray) and of the 
ellipsoids of type S2 (dark gray). The dashed curve contained in region 25n is the zero-pressure curve, while 
the one crossing regions 25n and 25i is the border of the existence region of the ellipsoids S3; the dashed 
curve in (b) is the border of the existence region of the S2. The Riemann ellipsoids of types I, II and III in 
the unshaded subregions of their existence regions are not elliptic and hence are Lyapunov unstable. The 
counterparallel S2 ellipsoids, which lie in the unshaded subregion of 25g 2 marked LS, are Lyapunov stable; 
the coparallel S2 ellipsoids in 25g 2 n Int(25m) are Lyapunov unstable (in figure (b), this part of region 23g 2 is 
'hidden' below 25m). 

Proposition 5 (i) All S-^-ellipsoids and all counterparallel ^-ellipsoids are Lyapunov stable. 

(ii) The coparallel ^-ellipsoids for b € 23s 2 n Int(Bni) (i.e., G(b±, 625^3) > 0) are not elliptic 
and hence are Lyapunov unstable. 

The first statement is due to Riemann |^(| , who proved it using arguments essentially equivalent 
to using the reduced Hamiltonian as a Lyapunov function. Riemann also showed that the reduced 
Hamiltonian has a saddle point, and thus cannot be used as a Lyapunov function, at any other 
ellipsoid. Riemann interpreted this fact as an indication that these ellipsoids are unstable ("labil"). 
However, as best we know this fact is unproven, except in the case of the S2 ellipsoids considered in 
statement (ii), the instability of which follows from Chandrasekhar's work, since he showed that the 
linearization of the equations of motion has an eigenvalue with nonzero real part at these points; 
see section 4.B. 

Remark: Given the structure of the problem, the Lyapunov stability of the equilibria of the 
reduced system corresponds to orbital stability of the solutions of the Dirichlet problem (see e.g. 

@, a ii)- 

The ellipticity of the Riemann ellipsoids has been studied by Chandrasekhar; most of his results 
are based on numeric calculations. We have repeated the numerical analysis of the ellipticity, 
finding the same results as Chandrasekhar for the S2 ellipsoids and essentially the same results for 
the ellipsoids of type I, but significantly different results for the ellipsoids of types II and III. The 
conclusions of our investigation are the following: 
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Numerical Conclusions 1: (i) The coparallel ^-ellipsoids for b in the complement of ¥>s 2 n 
Int('Bni), i-e. b satisfying 0(61,62,63) < 0, are elliptic. 

(ii) There are nonempty subregions of 23i, ¥>n and 2$m consisting of elliptic ellipsoids; these 
regions are the shaded regions in figures [| and 0. The ellipsoids of types I, II and III in the 
complement of these subregions are not elliptic and hence are Lyapunov unstable. 

This analysis has been carried out by analytically constructing the linearization of the reduced 
Hamiltonian vector field at the reduced equilibria and by numerically evaluating this matrix at 
some tens of thousands of ellipsoids of each type. Details about this analysis are given in the next 
subsection. (As explained there, Riemann sketches an argument that would imply statement (i).) 
Here we comment on these results. 

• The elliptic S2 ellipsoids of unknown Lyapunov stability lie between the "irrotational" curve 
and the border of region 23m • This is the dark shaded region in figure ||.b. 

• Our computations confirm the existence of two ellipticity regions for the ellipsoids of type I, 
as shown in figure 3 of 0, but disclose a finer structure of the lower region than that found 
by Chandrasekhar: the upper border of this region appears to transversally intersect the border 
b\ = 262 — 63 of the existence region at about 62/61 = .515, but there is a very narrow crescent 
of nonelliptic ellipsoids along this line, approximately in the range .5 < 62/61 < .503, which is 
shown in figure |4|.a. We could not determine with certainty if this crescent actually touches the 
line 61 = 262 — 63. 

• The region of elliptic ellipsoids of type II consists of three narrow fringes, shown enlarged in 
figure ||.b,c,d. None of these regions were detected by Chandrasekhar, who stated that no ellipsoid 
of type II is elliptic (0, pg. 169). Note that two of these fringes belong (for the most part) to the 
region of positive relative pressure considered by Chandrasekhar. 

• The region of ellipticity of the ellipsoids of type III has a three-lobed structure, the three 
lobes being separated by two very thin gaps; see figure ||.e, f. Chandrasekhar detected only the 
upper of these three lobes. To within the precision of our numerical computations, the two gaps 
close at approximately 62/61 = .01, but it is possible that, in fact, they persist to the origin. 

In an attempt to understand the discrepancies between our results and Chandrasekhar's, we 
tested our calculations using his linearization of the equations of motion, finding complete agree- 
ment. It is therefore possible that Chandrasekhar's conclusions were simply based on too few 
sample cases to obtain a detailed global picture of the structure of the regions of ellipticity. A more 
detailed comparison with Chandrasekhar's result is given in Appendix B. 

We conclude this analysis by observing that the sources of the richness and intricacy of the 
structure of the ellipticity regions resulting from these numerical computations are not yet un- 
derstood. Of course, we cannot be sure that our numerical computations, even though accurate, 
ultimately detected the finest details of this structure, particularly in view of the nontrivial struc- 
ture of the ellipticity regions for ellipsoids approaching a degenerate (flat or rodlike) configuration. 
Many features remain to be understood. For instance, figures ||.a and ^|f suggest that there are 
(nontrivial) relations between the ellipticity regions of the ellipsoids I (resp. Ill) and the boundary 
of the existence region of the ellipsoids S2 (resp. S3), which could be the object of a bifurcation 
analysis. Additional analytical insight into the problem is certainly needed. 

Remarks: (i) In his articles, Chandrasekhar refers to ellipticity as "stability" , while Riemann 
designated by this name what is now called Lyapunov stability. Because of this, and the previ- 
ously mentioned fact that Riemann regarded as unstable the ellipsoids which are saddle points of 
the Lyapunov function, Riemann and Chandrasekhar found different "stability" and "instability" 
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Figure 4: Details of the regions of ellipticity (shaded) of the ellipsoids of types I, II and III. All curves are 
as in figure [| 



Fassd and Lewis: Riemann Ellipsoids 



16 



regions. Chandrasekhar interpreted this difference as being due to some errors in Riemann's work 
(see page 875 in []|] and pages 185-7 in Q). Riemann's results are fully confirmed by our numerical 
computations (with the exception of the characterization of the saddle points as unstable, of course) 
and are, in fact, compatible with Chandrasekhar 's results once the distinction between Lyapunov 
stability and ellipticity is taken into account. 



(ii) In the quoted references, Chandrasekhar refers to [12| for an explanation of "the origin of his 
[Riemann's] errors" . In that article, to the best of our understanding, the source of these supposed 
errors is ascribed to the fact that Riemann treated the Dirichlet problem as a Lagrangian system 



(which is correct, as we have noted in Section 2), while according to [12] the system is Lagrangian 
only in the case of the ellipsoids of type S; remark (vi) in page 185 of || seems to indicate that 
Chandrasekhar shared this opinion. 

B. Numerical Computations. We now briefly describe the numerical procedures used in the 
verification of Numerical Conclusions 1 and in the Nekhoroshev analysis. We also provide the key 
components of the proof of statement (ii) of Proposition 5. We use canonical coordinates obtained 
by introducing a pair of local "Poincare coordinates" (q,p) on each sphere in the reduced phase 
space; these coordinates coincide with the Poincare elements used in the Kepler problem. The 
use of coordinates is not strictly necessary for the ellipticity analysis, but is very useful for the 
construction of the normal forms used in the analysis of the Nekhoroshev stability. 



Lemma 1 For any p > 0, the map 

(/ q 2 + p 2 / q 2 + p 2 q 2 + p 2 \ 

p y p — i — >~ q \l p — 4 — ' p — 2 — J ^ 4,1 ' ) 

is a diffeomorphism from the disk {{q,p) G M 2 : q 2 + p 2 < 2p} onto the half-sphere {to G M 3 : 
||to|| = p, 1713 > 0}. It is symplectic in the sense that (rh p )*(dp Adq)(m x u>,m x u/) = — m ■ uj x uj' 
for all to in the half sphere and all uj,uj' G M 3 . 

Proof. The map ( |4.1| ) is smooth and onto; its inverse to i— > y^2/(p + m3)(-m2,rai) is also smooth. 
If v = (vi,V2,vz) is tangent to the sphere at the point m = m p (q,p), then v = + 
for some constants v q and v p ; specifically, if we set k\ = 2p — q 2 — p 2 and &2 = k\ + 2/3, then 
kiy/k~2V q = qpv\ — (ki— p 2 )v2 and kwfk^Vp = {k2—q 2 )v\ — qpv2- Thus (rh p )^(dpAdq)(mxuj,mxio') = 
(m x uj)p (to x uj') q — (to x uj) q (to x uj') p = —to • {uj x lo') for any uj,lo' G K 3 . □ 

When studying a given equilibrium (6*,0, M*), with M* = (to ; *,to*), we shall translate the 
origin of the coordinates b to the equilibrium value b*; for simplicity, we will retain the notation b 
for the translated coordinates. Moreover, we shall use 'rotated' Poincare coordinates (qi,pi) and 
(q2,P2) centered at m\ and m* on the spheres; specifically, for given to* G M 3 we define 

(q,p) i ^ m(q,p;m*) := R m * rh\\ m *\\(q,p) 

where R m * is an orthogonal matrix such that i? m *e3 = Tj^snr- Since rotations are symplectic, 
the coordinates to(-;to*) are symplectic. These coordinates are not uniquely defined: right- 
multiplication of Rm* by a rotation about e% produces a new system of rotated Poincare coordinates 
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centered at m* . To compress the notation, we write (q,p) for {qi,q2,Pi,P2) andf] 

M(q,p;M*) = (rh(qi,pi; m*), m(q 2 ,p2; m*)) . 
Thus, the reduced Hamiltonian fl2.8| ) takes the form 

M(b,c,q,p;b*,M*) = ^c-K(b* + b)c + ^ M(q,p; M*) ■ J(6* + b)M(q,p; M*) + V(6* + 6) . 

We now study the block structure of the Hessian matrix !K"(0; b* , M*). To obtain a Hessian 
with as few nontrivial blocks as possible, we construct the coordinates using matrices such that 
Rm*£i = Rm*ei is parallel to a principal axis orthogonal to span[m;, m r ]. Using this choice of 
coordinates and denoting by "K bb the two-by-two block £ , etc, we have 



Lemma 2 For all equilibria, the blocks 3i" b , "K'cq, 3i'c P , an ^ •^qp e Q ua ^ zero an d ^'cc = 

Moreover, 'H'S = for the equilibria of type S. 

Proof. The equalities l H'{, b = "K" q = !H" p = and CH£' C = IK are obvious. Note that the 'standard' 



Poincare parameterization (44) of the sphere satisfies 



f)rn f)rh cP'rri 

— (0,0; p) = -^pe 2 , ^(0,0; p) = ^e u J^(0,0; p) = 0. (4.2) 

Using the latter equality and denoting by ^7 (resp. Jjp) the vector (resp. matrix) with entries 
given by the derivatives with respect to qj of the entries of M (resp. J), etc., we have 

d 2l K dM <91J 

_(0; 6 .,M.) = i i ^-(0,0).^( i .)M. 

7 f^-(0;6*,M*) = i? (0,0) -1(6*)^ — (0,0) (4.3) 
oqidpj dqi dpj 

d 2 "K dM <9J 

-(0;b*,M*) = R-^—(0, 0) • tt-Q)*) M* , 



dbidqj dqj db{ 

where -R is the six— by— six block diagonal matrix with diagonal blocks Rm* and Rm* • 

Due to the block structure of J, the two 'components' of J^-(6*)M* lie in span[m;, m r ]. (Here 
we refer to mi and m r as the 'components' of (mi,m r ).) On the other hand, ( [4.2j ) and the condition 
that R m *ei = R m *ei be orthogonal to span[m/,m r ] imply that the components of i?^-(0, 0) are 
orthogonal to span[m^, m r ] for j = 1,2; hence ( fOf ) implies that ^ ^ (0; b* , M*) = 0. To see that 
dq dp (fo ^*) e Q ua l s zero, note that R m *ei = R m *ei is, by construction, an eigenvector of Ji(b*) 
and J2(b*); hence ( |4.2| ) implies that the components of R T I{b*)R ^-(0, 0) are orthogonal to the 

components of §^-(0,0). 

In the case of the S ellipsoids, the two components of M* are both parallel to the same coordinate 
axis. Due to the block structure of J, the components of ^(b*)M* are parallel to the components 

5 In what follows, to keep the notation simple, we will write all vectors as row vectors, even though they must be 
regarded as column vectors in all matrix formulas. 
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of M*; thus, both components of R T -^-M* are parallel to e^. By ( |4.2| ) and (f4.3|), this implies 

•K, □ 

Let us now consider the ellipsoids of type S2 ■ With a suitable choice of the coordinates, the 
Hessian has the diagonal block structure 

X"(0;b*,M*) = diag^K,^,^]. 

We limit ourselves to stating without proof the following facts, which can be proven using elemen- 
tary, though sometimes laborious, computations: 

• The blocks K and "K qq are positive definite for all the S2-ellipsoids. 

• The block "Kp p is positive definite for all the counterparallel S2~ellipsoids and is indefinite for 
all the coparallel S2-ellipsoid. 

The block "K 1 ^ has a relatively simple expression, but the analytical study of its definiteness is 
quite laborious. We verified numerically that it is positive definite at all ellipsoids of type S2 • 
(Chandrasekhar's study of this block in Q is also numerical.) The linearization at the equilibrium 
of the Hamiltonian vector field of "K is X(b*,M*)(b,c,q,p), where the matrix X(b*,M*) is the 
product of the symplectic matrix and Di"(0; b* , M*). Thus, at a type S2 ellipsoid, the matrix 
X(b*,M*) is block diagonal, with the two 4x4 blocks 

Xbc = {-k b 0) and x ™ = ^ 

Since the eigenvalues of each of these two matrices are purely imaginary if its two nonzero blocks 
are positive semi-definite, we arrive to statement (i) of Numerical Conclusions 1. This conclusion is 
not rigorous because it is based on the numerical verification of the definiteness of !K^ fe . Statement 
(ii) of Proposition |5| follows from the definiteness of [Kg q and the indefiniteness of < K l ^ p . 

Remark: In [^] , Riemann suggests an analytic argument that would show that for all ellipsoids 
of type S satisfying to/ 7^ m r , i.e. &2> ^3) > 0, the function b 1— ► 5£(6, 0, 0, 0; b* , M*) has a 

minimum at b = 0. Since the block "K 1 ^ is positive semi-definite at such a minimum, this result 
could be used to construct an analytic proof of statement (i) of Numerical Conclusions 1. However, 



as Riemann explicitly states, he does not actually carry out the minimality argument in [20]. 
We now consider the ellipsoids of types I, II and III. The Hessian has the form 



K"(0;b*,M*) 



The eigenvalue problem of this matrix factorizes into the eigenvalue problems for the three matrices 

/ o/ff nrll 
nrll I ^bb ^bg 



( K'bb 





nrll 
K bq 







K 








n_r/l 
K qb 





nrll 
11 





V 








nrll 



One can prove that the two-by-two matrix "Kp p is indefinite for any ellipsoid of type I, II or 
III. (This implies that these ellipsoids are saddle points of the reduced Hamiltonian, as stated by 
Riemann.) The eigenvalue problem for the linearization of the Hamiltonian vector field determined 
by this Hessian does not factorize. Therefore the analysis of the ellipticity is particularly difficult. 
We know of no rigorous ellipticity results for these ellipsoids. 
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5 Nekhoroshev— stability of the Riemann— ellipsoids: Generalities 

We base our study of the Nekhoroshev-stability of the Riemann ellipsoids on some recent results 
in Nekhoroshev theory, which we describe in this section with reference to the problem at hand; 
complete proofs can be found in @ (see also @, 0, Q). As is common in Hamiltonian 
perturbation theory, the procedure consists of the construction of the Birkhoff normal form of 
order four for the reduced Hamiltonian and testing of the normal form for certain properties that 
we now specify. 



A. Birkhoff Normal Form. Consider an elliptic equilibrium (b*,M*) with M* = (rnf,m*), 
to* / / to*. Use the canonical coordinates £ = (b,c,q,p) defined as in Section 4.B. Let 
, M*) denote the j— th term in the Taylor series at £ = for the reduced Hamiltonian 
'K, so that "Kj is a homogeneous polynomial of order j in £ and 

M(£;b*,M*) =K 2 (£;b*,M*) + K 3 (Z;b*,M*) + W 4 (£;b*,M*) + .... 

Let X(b* , M*)£ be the Hamiltonian vector field of "K2- Since the equilibrium is elliptic, the matrix 
X(b*,M*) has purely imaginary eigenvalues 

±ioj j {b\M*), j = l,...,4, 

with the convention Uj > 0. 

The first step in the construction of the Birkhoff normal form for the Hamiltonian "K is the 
"symplectic diagonalization" of the quadratic term "K 2 - It is known that if the eigenvalues of 
X(b*, M*) are all distinct, then there is a linear canonical change of coordinates 

£^ S = [B, C, Q, P) = T(b*,M*)-^ 

such that 3< 2 (Z;b*,M*) := % 2 {T(b*,M*)E;b*,M*) has the form 

~ Bf + Cf bi + ci Qj + P? Ql + Pi 

5t2 = sxui h s 2 u\ h s 3 cj 3 h , 

where Sj = ±1. The matrix T and the numbers Sj are constructed as follows: For j = 1, ..,4, let 
xf = x'j±ix?! £ C 8 , with x'j, x] £ K 8 , denote any eigenvectors of the matrix X(b*,M*) associated 
to the eigenvalues ±icoj. One can show that 



Tj(b*,M*) := x'j-Jsx" + 0, j = l 



5 • • > ^ > 



where Jg = diag [J4, J4]; here J4 denotes the standard four-by-four symplectic matrix f ^ J V 
We set (uj,Vj) = (x'^x") if Tj > and (uj,Vj) = (x",x'j) if Tj < 0, and define 

/-iil «2 «1 V2 ^3 W 4 U 3 V A \ 

T = , Sj = SignH^), 

V71 72 71 72 73 74 73 74' 

where 7j = yjrTf; the notation means that the eight vectors are the rows of the matrix. Note that 
the condition that all the eigenvalues of X be distinct (and hence nonzero) can be regarded as the 
condition that the "frequency vector" 



= (s\UJl, . . . , S4W4) 
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has no resonances of order one or two: 

Q ■ v ^ for all v G I? such that \v\ = 1, 2, (5.1) 

where, for integer vectors, \i>\ = \v\\ + • • • + \ v$\. 

Remark: In the case of the ellipsoids of type S, the eigenvalue problem for the matrix X(b* , M*) 
factorizes into that for two 4x4 blocks. Since the construction of the matrix T also factorizes, it 
is necessary to check the nonresonance condition ([O]) for all resonances of order one, but only for 
the four resonances (1, ±1, 0, 0), (0, 0, 1, ±1) of order two. 

For the construction of the Birkhoff normal forms we use the complex coordinates U = 
(W 1} W 2 , W 3 , W 4 , Z u Z 2 , Z 3 , Z A ) defined by 

w iB i + c i ~ Bj + iCj iQ d + Pj Qj + iPj 

w > = -^^> Z > = ^VT^> W ^ = -^VT^> Z ^ = ^VT^ (J = 1 ' 2) 

(note that W are the coordinates, Z the momenta). Since 3 = T,U, with 

/ % 1 \ 
s 4 o\ . 1 % 1 

I arm > ,a = 

-1 -i 
\ -1 -ij 



the Hamiltonian in terms of the complex coordinates U is H(U;b* , M*) := ^K(TT,U;b* , M*) and 
its quadratic part is 

4 

H 2 {U;b*,M*) = ^2in j (b*,M*)Z j W j . 
j=i 

We now describe the construction of the Birkhoff normal form of order four of the Hamilto- 
nian H. Given a function / of W = (Wi, W2, W3, W4) and Z = (Z\, Z 2 , Z3, Z4), we denote its 
Taylor series by f(W,Z) = E^eN 4 fjkW j Z k , where W j = W^W^W^W^ etc. For any integer 
vector 1/ £ Z 4 , we define the f-th harmonic of / by 

(f) v (W,Z) = J] f jk WiZ k . 

j—k—v 

The average of / is its harmonic (/)o- The spectrum of / is Sp(/) = {v G : (/)„ 7^ 0}. 

We construct the Birkhoff normal form by means of the so-called Lie method. Let $ x denote 
the time-one-map of the flow of the Hamiltonian vector field of a function x\ the Lie transform 
generated by x °f a function / is / o If / and x are analytic in an open neighbourhood of a 
point (the equilibrium point, in our case), then / o is also analytic in some sufficiently small, 
but nonempty, open neighbourhood of such a point and there one has / o = Yl'jLo J! ^xf ■ (We 
adopt the sign convention Ljg = {f,g} = §^-§^ — ^~§w f° r the Poisson brackets.) 

The generating function x °f the Lie transform is constructed by solving the so-called homo- 
logical equation, which has the form {H 2 , x} = f ~ (f)o f° r some function /. The solution of this 
equation is formally given by x = where 

^es P (/)\{0} 
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Obviously, this is well defined if O does not resonate with any v £ Sp(/), and if the series converges. 
In our case / will always be a polynomial, either / = H 3 or / = H'^ (see below), so no convergence 
problems exist. 

The fourth-order Birkhoff normal form of H = H2 + H 3 + H± + . . ., if it exists, is constructed 
using two Lie transforms, which average the terms of degree three and four of H, respectively. The 
first Lie transform is generated by the solution 

Xi = $n(H 3 ) (5.2) 

of the homological equation {H2, Xi} = ^3 an d is well defined if 

Q(b*,M*) -u^O W £ Sp(H 3 ) . (5.3) 

Since (H 3 )q = and Jj^Hi = —L Xl H 3 one sees that this Lie transform conjugates H = H2 + H 3 + 
Hi + .... to 

H' = H 2 + H' A + ... , where H f A = ~L X1 H 3 + H 4 . 

The second Lie transform is generated by the solution \2 = §^(-^4) of {H21X2} = H'^ — {H' 4 )q. It 
is well-defined if 

n(b*,M*) -v^O W G Sp(#4) \ {0} (5.4) 
and leads to the fourth-order Birkhoff normal form 

H" = H 2 + H'l + ... with H'l = (±L X1 H 3 + H A ) Q . (5.5) 

For our purposes, the only important points are controlling the existence of the fourth-order 
Birkhoff normal form and testing that it possesses the convexity properties discussed above. Nei- 
ther of these tasks require the actual computation of the generator \2'- one needs only test ( j5"33| ) 
and compute H±. 

B. Nekhoroshev stability of elliptic equilibria. Let us introduce the action functions 
/ = (!]_,..., I4), Ij = iWjZj, and regard H'l as a quadratic form on M 4 by writing 

H'i{I) = h-A{b*,M*)I 

where A(b*, M*) is a 4 x 4 symmetric matrix. The Birkhoff normal form of order four (|5.5| ) is said 
to be 

• convex if the quadratic form H'[ is definite, i.e. the eigenvalues of A are either all positive or 
all negative. 

• quasi-convex if the restriction of H'l to the subspace orthogonal to f2 is definite, i.e. 

n ■ I = and H'l (I) = 1 = 0. 

• directionally quasi-convex if the restriction of the quadratic form H'l to the plane orthogonal 
to Q, is nonvanishing in the "first 16-ant" , i.e. 

n-I = 0, H'l(I) = and /i,...,/ 4 >0 =^ 1 = 0. 
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Clearly, each notion generalizes the previous one. 

If the Hamiltonian H(W, Z) and its Birkhoff normal form of order four are analytic, then any 
of these conditions implies the Nekhoroshev-stability of the equilibrium, in the precise sense that 
for any small e > one has 

||J(0)||<€ ||/(t)||<e Q for |t|<expe - 

for some positive constants a and (3. There are different possible values for the constants a and /?, 
which we report here with reference to the case under consideration of a system with four degrees 
of freedom: 

- If the Birkhoff normal form is directionally quasi-convex, then one can take a = j3 = \ and 
it is also possible to prove a stricter confinement of motions on correspondingly shorter times, 
namely a = |±|, /3 = ^ for any k > (see @). 

- If the Birkhoff normal form is quasi-convex (directional quasi-convexity does not seem to be 
sufficient), then it is possible to prove the above estimates with a = 1 and = jg. (See [pd|].) 

- If the Birkhoff normal form of order four is quasi-convex and it is possible to construct the 
Birkhoff normal form of order s > 4, then one also has a = 1, j3 = (see [pd]]). 



Remark: The original motivation for the quoted articles JTH], [11| on the Nekhoroshev stability 



of elliptic equilibria was precisely the present study of the Riemann ellipsoids. Until a few years 
ago, all Nekhoroshev stability results assumed that the frequency vector would satisfy a strong 
nonresonance condition, typically a diophantine one. This would be insufficient for the study of 
a problem such as the present one, in which the frequencies of the equilibrium, and therefore its 
nonresonance properties, depend on continuously varying parameters. The notion of directional 
quasi-convexity was introduced in [|| in connection with the study of the Nekhoroshev-stability 
of the triangular points of the restricted three-body problem. As will be seen below, this notion 
also plays a central role in the study of the Riemann ellipsoids, since the majority of the ellipsoids 
are not quasi-convex. (Nekhoroshev estimates also hold if a sufficiently high order Birkhoff normal 
form exists and satisfies some "steepness" condition; directional quasi-convexity allows us to avoid 
these generalizations, the hypotheses of which would be rather difficult to verify and which would 
lead to significantly worse values for the constants a and /?.) 

C. KAM theory. If the Birkhoff normal form of order four is nondegenerate, in the sense that 
det A(b* ,M*) ^ 0, then KAM theory applies, ensuring that in any sufficiently small neighbourhood 
of the equilibrium, the majority of the initial data gives rise to motions which are quasi-periodic 
with four frequencies. (See e.g. jjj Eq].) 



6 Nekhoroshev-stability of the Riemann ellipsoids: Numerical re- 
sults 

A. Results. We have numerically constructed the Birkhoff normal forms for a (quite large) 
number of sample Riemann ellipsoids. Specifically, within the existence region of each type of 
ellipsoid, we considered a mesh in the plane (j^> j^) determined by vertical lines uniformly spaced 
at a distance .0025; the number of mesh points on each vertical lines typically varied between 
twenty and one hundred, depending on the length of the line (however, in some cases we used as 
many as five hundred points). This analysis leads to the following 
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Figure 5: (a) The ellipsoids of type I that are resonant up to order four and (b) the quasi-convex Riemann 
Ellipsoids of type I in the lower region of ellipticity. 



Numerical Conclusions 2. The Birkhoff normal form of order four exists and is directionally 
quasi-convex and nondegenerate for all the Riemann ellipsoids except for those lying on a finite 
number of curves in the set ¥> corresponding to resonances of order up to four. 

We regard the fact that all of the computed nonresonant Birkhoff normal forms are directionally 
quasi-convex as indicating that all nonresonant Riemann ellipsoids are Nekhoroshev stable. We 
investigated the existence of low-order resonances in greater detail, without restriction to the points 
of the considered mesh (see subsection 6.B for some details), finding that there are a finite number 
of curves in the plane (j^-, at which at least one of the nonresonance conditions ( |5,1|) , ([5.3j), 



(5^) is violated. Specifically, we found 8 different resonances for the type S2, 52 for the type I, 33 
for the type II, and 47 for the type III. For example, figure [|.a reports the resonant curves for the 
ellipsoids of type I. 

We found that the Birkhoff normal form is quasi-convex only for a few ellipsoids of each type; 
in turn, very few of these ellipsoids are actually convex. For example, figure ||.b shows the quasi- 
convex ellipsoids in the lower region of elliptic type I ellipsoids. The determination of the quasi- 
convex ellipsoids has some interest because, as discussed in the previous Section, such equilibria 
may have stronger stability properties if they are nonresonant up to sufficiently high order. (From 
this perspective, the finer distinction between convexity and quasi-convexity does not instead seem 
to be as significant.) 

Nondegeneracy of the Birkhoff normal forms implies that the KAM theorem applies, so that 
the majority of the motions near each equilibrium in the reduced phase space are quasi-periodic 
with four frequencies. "Reconstruction" of these motions then gives quasi-periodic motions of the 
Dirichlet problem, with up to (and in fact, typically) eight frequencies. 



B. Numerical Procedures. We give now some information about the numerical procedures 
adopted here. The fourth-order Taylor expansion of the Hamiltonian in coordinates (b,c,q,p), 
and hence the Hamiltonian matrix X(b* , M*), has been constructed analytically for each type of 
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ellipsoid. The coefficients of these polynomials have been evaluated for any ellipsoid in our mesh, 
and all other operations have then been performed on polynomials with numeric coefficients. These 
include: 

- Computing the eigenvalues and the eigenvectors of X(b*,M*) and testing for the presence 
of resonances of order one and two. Note that once the eigenvectors are known, the matrix 
S(b*,M*) determining the coordinates H and U is also known. 

- Constructing the normal forms by means of formulas ( |5.5D and ( |5.2| ) . This requires computing 
Poisson brackets, i.e. derivatives, and averages of polynomials with numeric coefficients; these 
are easily implemented operations. 



- Checking the nonresonance conditions fl5,3|) and (5.4), which assure that the normal forms can 
be constructed. This involves determining the spectrum of a polynomial. 

- Testing the normal forms, namely the numeric vector £1 and matrix A, for convexity, quasi- 
convexity, or directional quasi-convexity. 

The chosen implementation of the latter two points requires some explanation. 

Spectra and Resonances. We computed the spectra of H% and H'l at every point (&*, M* ) in the 
considered mesh and took the unions S3 and £4 of all these sets. (The union is taken because the 
numerically computed spectrum at a given point may lack very small harmonics.) We identified as 
resonances those v G £3 U X4 for which the scalar product $7(6*, M*) ■ v either vanishes or changes 
sign on our grid; of course, use of a finer grid could reveal additional resonances. (The resonant 
curves of figure ||.a were constructed by computing the approximate zeros of 0(6*, M*) ■ testing 
for changes of sign on a finer grid.) 

Quasi and Directional Quasi- convexity. If the matrix A is not convex, then the tests for quasi- 
convexity and directional quasi-convexity can be carried out as follows. Let R £ S0(3) be such that 
R£l = (1,0,0,0), so that the Hessian of the restriction of the quadratic form H'l to the subspace 
orthogonal to O is the lower right hand 3x3 block A of the matrix RAR T . Quasi-convexity of H'l 
is equivalent to convexity of the matrix A. 

Under the additional hypothesis that the restriction of H'l (I) to the subspace orthogonal to 
is nondegenerate|] the test for directional quasi-convexity can be performed as follows. The 
quadratic form H'l is directionally quasi-convex iff any "asymptotic vector", namely any nonzero 
vector / satisfying / • AI = 0, points out of the first 16-ant, i.e. has at least one negative and 
one positive entry. Since these vectors form a cone (of dimension two, if A is nondegenerate) we 
need only test the asymptotic vectors in the intersection of the cone and the unit sphere, namely 
two ellipses. In order to determine these ellipses, note that in the absence of quasi-convexity, the 
matrix A has at least one eigenvalue of either sign. We may assume that the eigenvalues of A 
satisfy at > and 02,03 < 0, replacing A with —A if necessary. (This is possible because A 
and —A have the same asymptotic vectors. Also, the nondegeneracy hypothesis assures that all 
of the eigenvalues are nonzero.) Let {^1,^2,^3} be an orthonormal eigenbasis associated to A, so 
that S = (v\,V2, V3) € SO(3) satisfies S T AS = diag ( | «i | , — 1 02 1 , — | «3 1 ) • In this basis, the cone of 
asymptotic vectors is given by 

(x = (xi, X2,Xs) £ I 3 : I? = -^J-x^ + ^~~^ a:; i)' • 
6 This hypothesis is not crucial, but it is satisfied by all the computed ellipsoids. 



Fassd and Lewis: Riemann Ellipsoids 



25 



The two ellipses of unit vectors on this cone consist of the vectors 

x^iO) = \±J ^ , cos 2 6> + ^ , sin 2 fl, ,/ a \ . cose, J a \ . sing) , 9 £ S 1 . 
\ yai + |a2| ai + |a 3 | y a.\ + |a 2 | y «i + |ck 3 | J 

In the original basis in the subspace orthogonal to S7, these vectors are 

i ± (9) = (i 1 ,i 2 ,i 3 ) = sx^e). 

These are three-dimensional vectors belonging to the subspace orthogonal to (1,0,0,0) = RQ, 
so the corresponding four-dimensional vectors are 1^(9) = R T (0, I\, I2, 13). The normal form is 
directionally quasi-convex iff for every 9 € S 1 the vector I + {9) has at least one negative and 
one positive entry. (It is obviously sufficient to consider only one of the two circles.) In the 
computations, we have tested this condition by checking that, at each of a number (typically 101) 
of uniformly spaced values of 9 in the interval < 9 < 2tt, there are two eigenvalues whose product 
is negative. 



7 Conclusions 

In this work, we have posed the problem of the Nekhoroshev stability of those Riemann ellipsoids 
that are elliptic but of unknown Lyapunov stability. We have provided numerical evidence that 
all Riemann ellipsoids that are not resonant up to order four are Nekhoroshev stable. We have 
developed a consistent and rigorous Hamiltonian formulation of the problem on a four-to-one 
covering manifold, which validates the procedures (and clarifies the presence of a previously not 
understood discrete symmetry in the problem). Within this formulation we provided explicit, 
precise descriptions of the existence regions and of the momenta of the Riemann ellipsoids of the 
various types and we reviewed the existent results on the stability and ellipticity of the Riemann 
ellipsoids, showing (numerically) that some of the latter regions are significantly larger than was 
previously thought. 

At the conclusion of this work, it is our opinion that several key features of the problem are not 
yet understood. Some open problems include: 

- Properties of axisymmetric Riemann ellipsoids other than the Maclaurin spheroids. 

- Analytic determination of the regions of ellipticity. Understanding of their subtle structure. 

- Lyapunov stability or instability of those Riemann ellipsoids that are saddle points of the 
reduced Hamiltonian. 

Of course, it would also be interesting to prove the Nekhoroshev-st ability rigorously, that is, without 
recourse to numerical calculations, but this appears to be a formidable task at present. 

In conclusion, we would like to mention a point of historic interest. Riemann's treatment of the 
Dirichlet problem is exceptionally advanced and uses several geometric constructions, including 
Poisson reduction, the amended potential, and the trivialization of the tangent and cotangent 
bundles of Lie groups, that are now considered to be central to modern geometric mechanics. 

Riemann's apparent lack of concern with proving the actual instability of the saddle points of 
the reduced Hamiltonian is probably typical of that era — this error permeated the mechanical 
literature long after Lyapunov clarified the subject at the beginning of the twentieth century. At 
the dawn of a new century, the Lyapunov (in) stability of equilibria that are both elliptic and saddles 
of the 'obvious' Lyapunov function is still a formidable problem. Nekhoroshev stability provides a 
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weaker, but very practical, notion of stability, which has the ultimate advantage of being compatible 
with Lyapunov instability. 



8 Appendix A: On the proof of Proposition g 

The proof of Proposition ||, even though elementary, is too lengthy to be reported here. Hence we 
indicate the main steps in the form of Lemmas that we do not prove. The relative equilibria are 
the critical points of the reduced Hamiltonian H„ : P v — > R and hence are given by the conditions 

V b H v = 0, V c H v = 0, mxV m H v = 0. (8.1) 

Since VcH^ = iff c = 0, we can equivalently search for critical points (b, m) of 

H v (b,m) := ±m-J(6)m + V(b) . 

Note that only the first equation in flS.ip depends on the potential V. 



Lemma 3 A nonzero vector m = (mi,m r ) 6 (R ) satisfies m x X7 m H v (b,m) = iff either of the 
following alternative conditions is verified: 

(i) mi and m r are both parallel to the same principal axis; 

(ii) mi and m r both belong to the same principal plane, say e« © ej with bi > bj, but are not 
colinear with either principal axis, and, together with b, they satisfy any of the following three 
conditions (where (i,j,k) is a permutatior^\ of (1,2, 3)): 

(ii. 1) bj < | bi — 2bk \ and 

ij - m\.\2 _ ( bj + 6 fc \4 b\ - (bj - 2b k ) 2 



m\ + m\.J \bj-b k J bf - (bj + 2b k ) 2 
m\ — m 3 r m] + m l r (bj + b k \ 2 (h + b k \ 2 h + bj — 2b k 



m\ + m! r m\-m\\ bj - b k 



/ \bi-bj h + bj + 2b k ' 1 ' ' 



(11.2) bi = 2b k + bj, m\ = —m l r , and mj 

(11.3) bi = 2b k — bj and mi = m r . 



m 



Equations ( |8.2|) and (|8.5| ) imply that ( |8.2[) is also valid with the index i and j exchanged. If 
either mi or m r vanish (the case of an irrotational ellipsoid), then equations ( [3,2; ) and (|8.3| ) are 
equivalent to k = 2 and b\b\ + b\b\ + b\b\ = 36|. 

We now consider the equation VbH ri = 0. Set 

% = {bG'B : G^(bi,bj,b k ) > , G% (h, bj, b k ) > 0} , j = 2,3 

where (j, k) is any permutation of (2, 3) and 

S± = {b e S : bj < ±(bi - 2b k ) , D(bi, bj, b k )^0, G*(b it bj, b k ) > , G^(bj,bi, b k ) > 0} 

where (i,j, k) is any permutation of (1, 2, 3) with i < j. (We shall see below that some of these sets 
are empty.) 



7 By a permutation, we mean either an even or an odd permutation. 
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Lemma 4 Assume m = (mi,m r ) £ (R 3 ) 2 is nonzero and satisfies mx'V rn H ri = 0. Then V&i^ = 
iff one of the following two conditions is satisfied: 

(i) b £ 'Bij, mi = rnjej, m r = miej and (mj ± mi) 2 = G±{b\,bj,bk) for either j = 2, k = 3 or 
j ~ 3, k - 2. 

(ii) b e r B~ j ] J'Bf j , mi = mlei+mjej, m r = mlei + miej, {m\±m % r ) 2 = G±(bi,bj,bk), (mj±mi) 2 = 
G^(bj,bi, b k ) and, if bi / 2b k ± b j; 

I , (mi + ml)(m) + m l S)\ i ~ ~ 

m - mi = sign (2b k -b % - b^ ¥±-L r ± JG^bj,^, b k ) , 

V m\ — m l r j v 

for any permutation k) of (1, 2, 3) u;zf/i i < j. 

The proof of part (i) of Proposition [| is completed by showing that the "existence regions" of 
the previous Lemma coincide with those of the Proposition and that all the equilibria described in 
that Lemma are Z4 x Z2-equivalent or adjoint to those indicated in the Proposition: 



Lemma 5 The sets B^j, T> 23 and ¥> 12 are empty, while ¥>2 = 25s 2 , S3 = 25s 3 , 25 13 = 25i, 25j" 3 = Sn, 
and 15+ = S m . 



Lemma 6 Under the assumptions of Lemma 0, a poini (6, 0, m) is a relative equilibrium for some 
reduced system iff it is Z4 x ^-equivalent to one of the points (b, 0, (//+ (6), /i+ (&)))■ 

Part (ii) of the Proposition is a consequence of the following 



Lemma 7 TTie maps /ij, a = S3, II, III, are never zero. For a = S2 and a = I, i/iey vanish on the 
sets indicated in Proposition 



9 Appendix B: On Chandrasekhar 's linearization 

To facilitate the comparison of our results with Chandrasekhar's, in figure || we display the regions 
of elliptic ellipsoids of types II and III using the coordinates employed by Chandrasekhar, namely 
(l^'ftf) f° r * ne ellipsoids of type II and (5^,^) for the ellipsoids of type III. (For the latter, 
Chandrasekhar also uses the coordinates the corresponding figure should be compared 

to figure ||.b.) Among all these ellipsoids, Chandrasekhar identified as elliptic ("stable", in his 
terminology) only the ellipsoids in the lowest ellipticity region in figure ||.d. 

As we have already remarked, Chandrasekhar's linearization of Riemann's equations yields 
the same results as ours (after correcting a few misprints that we report in the Remark below 
for the convenience of the interested reader). Hence, it seems probable that the meshes used by 
Chandrasekhar in his numerical study were too coarse, or too small, to capture all of the ellipticity 
regions shown in figure 0. Certainly, the very few ellipsoids whose frequencies are reported in H 
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b3 :b2 b 3 ±2 




Figure 6: Regions of elliptic (dark gray) and non-elliptic (light gray) ellipsoids of types II (a-c), shown at 
three different resolutions, and III (d). The indicated individual points are the points tabulated in [gj. The 
dashed curve in (a-c) is the zero pressure curve. 

(three of type II and four of type III, which are indicated by large points in figure ^) are consistent 
with our result. 

Chandrasekhar did not exploit (nor, apparently, notice) the Lagrangian and Hamiltonian struc- 
ture of the Dirichlet problem and worked with the full sixteen-dimensional Riemann's system, 
without exploiting the symmetry of the problem. The linearization of the reduced and unreduced 
systems are, of course, equivalent: the four additional frequencies appearing in the full system (two 
of which are identically zero) do not alter the ellipticity. 

Remark: In §, the frequencies of the linearized equations of motion are determined by the 
system of nine equations (116). The nine- by-nine determinant of the coefficients of this system is 
given in formula (133) of [||], where there are two misprints: the entry (7,2) should be multiplied by 
2 + ((3 — 2^a\ja\ ) and the entry (8,3) should be multiplied by 2 + (7 — 2)[5a\ ja\. This determinant 
is reduced to the eight-by-eight determinant (134), in which the entry (1,8) should be changed to 
n 2 n 3 + 2^n 2 n\, the entry (2,7) to n 2 n 3 + 2^n\n 3 and the entry (3,5) to 2 ^3 + ^l n h the 
entries (8,1), (7,2) and (5,3) should be changed accordingly: the entry (i, j) is obtained from the 
entry (j, i) by exchanging Qf. and fij. (k = 2,3). Also, the symbol -B123 appearing in equation (135) 
of ||, should be defined as in equation (104), Chapter 3, of p. Let us explicitly note that these 
misprints are not the sources of the discrepancies between Chandrasekhar's and our conclusions: it 
is only after they have been corrected that equation (134) of M detects as elliptic the regions that 
Chandrasekhar describes as "stable". 
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